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ABSTRACT 

The gravitational lens equation for the tilted Plummer family of models can be reduced 
to a one-dimensional non-linear equation. For certain values of the slope of the radial 
profile it can be reduced to a polynomial form. Both forms arc advantageous to find 
the roots, i.e. the images for a given model. The critical curve equations can also 
be reduced to a non-linear or polynomial form, and therefore it is useful to find the 
caustics. This lens model family has ample use in gravitational lens theory, and can 
produce up to five images. 
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1 INTRODUCTION 

Smooth non-singular isothermal sphere (SNIS) models are 
often used in gravitational lens theory. Including ellip- 
ticity, these models are call ed elliptical SNIS (ESNIS) 
Blandfor d fc Kochanekl Il987l ) or tilted Plummer family 



Kassiola fc Kov ncr 1993). This class of models can be gen- 
erated from elliptical mass distributions or elliptical po- 
tentials. The models obtained by means of mass distribu- 
tions are realistic even for higher ellipticities. The mod- 
els generated using the elliptical potentials are restricted 
to small ellipticities (Kassiol a fc Kovnerlll993l ). because the 
isodensity contours become dumbbell-shaped for higher el- 
lipticities. Moreover, the density for these elliptical poten- 
tials could be negative for cert ain combinations of parame- 
ters (|Kassiola fc Kovner|[l993l ). For more information about 
these models, the interested reader may consult the refer- 
ences at the end of this Letter. 

For ESNIS models, the gravitational lens equation (GLE) 
becomes non-linear, and hard to solve if the external per- 
turbations (density and shear) are taken into account. How- 
ever, for ESNIS models obtained from elliptical potentials, 
it is possible to simplify this equation to a one-dimensional 
non- linear or polynomial form. We will show how to do it. 
The advantage of reducing the GLE to a one-dimensional 
non-linear or polynomial form is that there is simple itera- 
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tive methods to get the roots of this class of equations. The 
critical curve equation can also be reduced to a non-linear or 
polynomial form. This allows us to find the caustics easily. 
Examples of these models will be presented in this Letter. 



2 ESNIS MODELS 

The tilted Plummer family of elliptic potentials has the form 



/ 2 i 2 i 2 



■lf-i> , < /3 < 1 

P = o, (1) 



where x = (xi, X2), x\ and £2 are measured in the direction 
of the principal axes of the ellipsoid (image position), j3 is 
the slope of the radial profile of the potential, xo is the core 
radius (same units as xi and X2), Ko is a constant, e, = lie 
(i = 1, 2) with e as the lens ellipticity, and ipQ is a constant 
potential reference, which will be taken as null. 
The scaled deflection angle components for the ESNIS mod- 
els have the following form: 



*(*) = 



(X% + £\X\ + 62 



xiy-v 



(2) 



where k% = «o£s- 

The surface mass density is given by 

k(x) = n (xl + six* + £2^2) l x o + PiSiX* + P2S2X2], 

where = 1 + {fi - l)ej. 

To avoid negative values of surface mass density and the 
dumbbell form of the isodensity contour shape, the ellipticity 
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e sho uld be restricted to small values (|Kassiola fc Kovnerl 
Il993h : 



3 



3-/3' 

Among the models, that the ESNIS class contains are: 

• power- law mass (e = 0, and xo = 0), 

• NIS (,3 = 1/2, andesC 1/5), 

• Plummer (/3 = e = 0). 



3 REDUCTION TO A NON-LINEAR OR 
POLYNOMIAL FORM 



The scaled ray tracing equation, or GLE is l|Schneider et al.l 
1 19921 ) 



y = Mx — a(x), 



(3) 



where y is the source position on the source plane, and the 
matrix M is given by 



M = 



1 — (j — 7 cos 20 —7 sin 20 
—7 sin 2(f) 1 — a + 7 cos 2<f> 



(4) 



where er, and 7 are the density and shear of external pertur- 
bations (nearby galaxies or a cluster contribution), and cj> is 
the shear angle. 

Let us define some variables for the sake of simplicity 



X = xl + E\x\ + £2X2, 



(•>) 



/(x) = x 1 -^x(/) = / 1/(1 -«. 

With these definitions the GLE system takes the form 



2/i 



U-2 



Mh - y I + M\2.'-'. 
M22 - y ) + l/_,j.r t- 



Solving for x\ and X2, we get 



t<2 



gxi = ( M22 - -j ) yi - M122/2 =: fi, 

gX2 = (Mil - y ) '/J - V 21 ;/i =: /_>. 
where 

</ - ( A/: 1 y) ( A/.-, y ) - 

Substituting © and © into (J5]) yields 
3 2 X(/) = ^off 2 + ei/i 2 + £2/2- 



(6) 
(7) 

(8) 
(9) 

(10) 

(11) 



By substitution of the right-hand side of equations ((8j), (|9]), 
and (|10[) into the last equation, we obtain 



(x(/)-^)Pi(/) = p 2 (/), 

where 



(12) 



Pi(f) 



and 



[fgf (13) 
[(Mu/ - K1XM22/ - fa) - Mi 2 M 21 / 2 ] 2 
[A-f/ 2 - (K1M22 - K 2 Afn)/ + kik 2 ] 2 , 



P 2 (/) = {ei[(M 22 /-K 2 )yi -M 12 j/ 2 /] 2 

+ e 2 [(M 11 /-K 1 )y 2 -M 21 j/ 1 /] 2 }/ 2 

= (A^-AxJ + Ao)/ 2 , (14) 

with 

.M = det[M] = MnM 22 - M 12 M 2 i, 

A 2 = ei(M 22 yi - M 12 y 2 ) 2 + ea(Muya - M 2 ij/i) 2 , 

Ai = 2[K 2 eiyi(M 22 j/i - M i2 y 2 ) + Kie 2 y 2 (Mnj/ 2 - Mnyi)] 

and 

A = £l(K 2 J/i) 2 + £ 2 (kiJ/ 2 ) 2 . 

Pl(/) and P 2 (,f) are both fourth order polynomials. 
From (fl2l we get 



/[A(/)] 



1-/3 



tft(/)+xoft(/)] 



1-/3 



(15) 



The last equation is a non-linear one-dimensional equation 
for / that can be solved by means of it erative numerical 
meth ods, for example, the Brent method (|More fc Cosnardl 
119801 ). 

Now, if we set j3 = m/n, with m, n integers (n ^ m), then 
(|15|) becomes a polynomial form: 



riPiifW 



[ft(/) + 4fl(/)] r 



(16) 



The image positions can be obtained as follows: first of all 
solve (|15p or (|16|) ; now substituting the /-values into the 
equations (H), JS), and (|10|) yield /1, / 2 , and 3; finally use 
equations (|5J and ([9]) again to get £1 and x 2 - 



4 THE CRITICAL CURVES AND CAUSTICS 

The critical curve is the curve formed by all image positions 
on which the determinant of the Jacobian vanishes and the 
caustic curve is the projection of this critical curve to the 
source plane. 

The components of the Jacobian J = dy/dx are 



J11 



dyi 
dx\ 

dyi 



Mn-j) +2(l-P)Kie 1 x 2 lX rl 



Jl2 = ^- = M12 + 2(1 - (3)K 1 £2X 1 X2X\ 
OX2 



J22 



dy2 

3X2 



y ) + 2(1 - 3)k2S2X 2 2X V 



where 77 = /3 — 2. 

The determinant of the Jacobian is 
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+ 

x 
+ 



J = det[J] = JuJ-a- JiaJai (18) 

= I7aW [((Mii/_Ki)Pi(/) 

2(l-/3) Kl£l / 3 [C 1 /-Ac 2 y 1 ] 2 x") 
((M 22 /- K2 )A(/) 

2(l-/3) K2£2 / 3 [C 2 /- Kl y 2 ] 2 X ") 

- (2(1 - /8)«ie a (Ci/ - K 2 yi)(C 2 f - «uy 2 )/V 

- M 12 /Px(/)) 

x (2(1 - /3) K2£l (Ci/ - « 22 /i)(C 2 / - ki»s)/V 

- M 21 /P!(/))] , 

where 

Ci = M 22 2/i - Mi 2 y 2 
and 

C 2 = Mii2/ 2 - M 2 iyi. 

In (fT8|) we have substituted ((Sj) and ©. From (fT8)> . it is clear 
that the critical and caustic curves could also be found by 
solving a one-dimensional non-linear equation or a polyno- 
mial in /. 

5 EXAMPLES 

5.1 Plummer's model 

Let us take (3 = 0. We get from (|15p the following polyno- 
mial: 

Bsf - P 4 / 4 + B A f + B 2 f 2 + Pi/ - Bo = 0, (19) 

where 

Ps = M 2 , 

B4, = M 2 xl + eid + e 2 Ct + 2C Z M, 

B 3 = 2k 1 k 2 M + {2x'qM + C 3 )C 3 + 2K 2 e 1 y 1 C 1 + 2«: 1 e 2 y 2 C 2 , 
P 2 = 2k x k 2 Cz - K? 2 E\y\ - K 2 e 2 y 2 - x\\C\ + 2n 1 K 2 M], 
Pi = KlK2[KlK2 — 2xqCs] 

and 

Po = [xokik 2 ] 2 , 
with 

C 3 = mM 22 +K2M11. 

The Source, images, critical curves, and caustics for this 
model are shown in Fig. 1. The values of the model param- 
eters are a = 0.05, 7 = 0.2, <j> = 36°, £ = 0, ko = 1 and 
xq = 0.5. The critical curves and caustics are inclined, be- 
cause we have chosen 0^0. 

5.2 ESNIS model 

Let us take (3 = 1/2. We also get from (|15p a polynomial: 
Bef - Ps/ 5 
where 
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Figure 1. Source (large black dot near the centre), two extended 
images (one of them is a long blue arc) and one tiny image at 
the centre (green dot), the critical curves (magenta ellipse and 
small circle), and the caustics (black diamond and dotted ellipse) 
produced by the Plummer model for following parameter values 
(3 = e = 0, 7 = 0.2, a = 0.05, <j> = 36°, k = 1, x = 1/2. 
The green dots are the positions of the images for a point source. 
The position of the extended source is (0.65, 0.45). The critical 
curves and caustics are inclined, because <j> ^ 0. The window size 
is 500 x 500 pixel. 



2C 3 M 2 , 

C| + 2k iK2 M - [M 2 x 2 + eid 2 + e 2 C|], 
2[Cz(x 2 ) M - «i« 2 ) + K 2 eiyiCi + Ki£ 2 y 2 C 2 ], 



k 2 £i2/? - K 2 e 2 y 2 , 



B 4 f 4 - P 3 / 3 + P 2 / 2 + Pi/ - Po = 0, (20) 



B 6 

Br, 
P 4 
P 3 

P 2 = [kik 2 ] 2 — Xq[2kiK 2 M + C, 

Pi = 2K\K 2 X%C'i 

and 

Po = [x KlK 2 ] 2 . 

In Fig. 2 the source, images, critical curves, and caustics for 
this model are shown. The values of the model parameters 
are a = 7 = 0.05, (j> = 36°, e = 0.2, = 1, and xq = 
0.5. The critical curves and caustics are again inclined ((f) 7^ 
0). For this model, one gets up to five images. We have 
tried with other /3-values (with different values for the other 
parameters) to find more than five images, but all of our 
ESNIS models always produce no more than five images. 
We have designed XFGLENSES, an interactive pro- 
gram intended to v isualize and model gravitational lenses 
ijFrutos-Alf aro 2001). The modelling part of this program is 
not finished yet. The Figures 1, and 2 were generated with 
this software. The interested reader can download it from 

\protect\vrule widthOpt\protect\href {http : //www. tat . physik.uni-tu( 

On this Website, there is also extensive information about 
the software. To solve the GLEs, the Brent method was im- 
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Figure 2. Source (large black dot near the centre), five images 
(one of them a long blue arc and one small image near the source), 
the critical curves (magenta ellipse and small circle) , and the caus- 
tics (black ellipse and diamond) produced by a ESNIS model for 
following parameter values f) = 1/2, e = 0.2, 7 = a = 0.05, (f> = 
36°, ko = 1, %0 = 1/2. The green dots are the positions of the 
images for a point source. The position of the extended source is 
(0.2, 0.2). The critical curves and caustics are inclined, because 
0^0. The window size is 500 X 500 pixel. 

plemented in the program. To visualize t he extended images, 
we use the Kayser-Schramm m ethod (jSchramm fe Kavserl 
ll987l ; lKavser fc Schrammll 19881 ). 



6 CONCLUSIONS 

Although the ESNIS model generated from elliptical poten- 
tials have an ellipticity restriction, these models are very 
useful in gravitational lens theory. The polynomial reduc- 
tion presented here is useful not only to invert the GLE, 
but also to find the critical curves and caustics of a given 
ESNIS model. Moreover it could be useful to model gravi- 
tational lenses with small ellipticities. This class of models 
always produces up to five images. XFGLENSES can be 
used to visualize the results of this class of models. To see 
the capabilities and models that this software has available, 
visit the Website mentioned above. 
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